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TOLLMIEN-SCHLICHTING/ VORTEX  INTERACTIONS  IN 
COMPRESSIBLE  BOUNDARY  LAYER  FLOWS 


Nicholas  D.  Blackaby1 
Department  of  Mathematics 
University  of  Manchester 
Manchester,  M13  9PL 
United  Kingdom 


ABSTRACT 

The  weakly  nonlinear  interaction  of  oblique  Tollmien-Schlichting  waves  and  longitudinal 
vortices  in  compressible,  high  Reynolds  number,  boundary-layer  flow  over  a  flat  plate  is  con¬ 
sidered  for  all  ranges  of  the  Mach  number.  The  interaction  equations  comprise  of  equations 
for  the  vortex  which  is  indirectly  forced  by  the  waves  via  a  boundary  condition,  whereas 
a  vortex  term  appears  in  the  amplitude  equation  for  the  wave  pressure.  The  downstream 
solution  properties  of  interaction  equations  are  found  to  depend  on  the  sign  of  an  interaction 
coefficient.  Compressibility  is  found  to  have  a  significant  effect  on  the  interaction  properties; 
principally  through  its  impact  on  the  waves  and  their  governing  mechanism,  the  triple-deck 
structure.  It  is  found  that,  in  general,  the  flow  quantities  will  grow  slowly  with  increasing 
downstream  co-ordinate;  i.e.  in  general,  solutions  do  not  terminate  in  abrupt,  finite-distance 
‘break-ups’. 


'This  research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Con¬ 
tract  No.  NAS  1-19480  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in 
Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23681. 
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1.  INTRODUCTION 


The  nonlinear  interaction  between  two  oblique  three-dimensional  Tollmien-Schlichting 
(TS)  waves  and  their  induced  streamwise  (longitudinal)- vortex  flow  is  considered  theoret¬ 
ically  for  a  compressible  boundary-layer  flow;  this  study  is  an  extension  of  the  paper  by 
Hall  &  Smith  (1989)  who  considered  an  incompressible  boundary-layer  flow.  The  same 
theory  applies  to  destabilisation  of  an  incident  vortex  motion  by  sub-harmonic  TS  waves, 
followed  by  interaction.  The  interaction  is  considered  for  all  ranges  of  the  Mach  number 
in  order  to  investigate  the  effect  of  flow-compressibility. 

The  motivation  for  such  a  study  is  essentially  the  same  as  expressed  by  Hall  &  Smith 
in  the  introduction  to  their  paper;  namely  that  often  in  experimental  studies  of  laminar- 
to-turbulent  transition  on  a  flat  plate  (eg.  Aihara  &  Koyaina,  1981;  Aihara  et  al ,  1985), 
there  appear  to  be  longitudinal  vortices  co-existing,  and  interacting,  with  the  viscous  TS 
modes.  As  there  is  no  concave  curvature  of  the  surface,  these  longitudinal  vortices  are 
not  Taylor-Gortler  vortices  driven  by  surface-curvature  (see  Hall,  1982a, b  and  subsequent 
studies);  instead  one  could  postulate  that  they  are  in  fact  being  driven  by,  and/or  inter¬ 
acting  with  TS  modes.  The  reader  is  referred  to  the  paper  by  Hall  &  Smith  (1989)  for  a 
fuller  account  of  relevant  experimental  findings,  as  well  as  supporting  computational  work 
(see,  for  example,  Spalart  &  Yang,  1986).  These  experimental  studies  are  all  for  incom¬ 
pressible  flow;  the  author  is  unaware  of  any  experimental  work  specifically  relevant  to  this 
compressible  study.  We  note  in  passing  that,  for  compressible  flow  over  a  heated  plate, 
buoyancy-driven  longitudinal  vortices  may  also  be  possible  (see  Hall  &  Morris,  1992). 

Recently,  the  origin  of  streamwise  vortices  in  a  turbulent  boundary  layer  has  been 
investigated  theoretically  by  Jang  et  al  (1986).  The  Reynolds  number  is  taken  to  be 
finite  and  their  formulation  is  of  the  Orr-Sommerfeld-type.  They  show  that  two  oblique 
travelling  waves  can  combine  nonlinearly  to  produce  a  stationary,  streamwise  vortex  — 
this  is  essentially  the  theoretical  idea  later  used  by  Hall  &  Smith  in  their  independent  work. 
However  the  latter’s  approach,  that  adopted  in  this  paper,  takes  advantage  of  the  feature 
that  the  Reynolds  numbers  of  interest  in  reality  are  large  and  so  the  Reynolds  number 
is  taken  as  a  large  parameter  throughout.  The  nonlinear  interaction  is  powerful,  starting 
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at  quite  low  amplitudes  with  a  triple-deck  structure  for  the  TS  waves  but  a  large-scale 
structure  for  the  induced  vortex,  after  winch  strong  nonlinear  amplification  can  occur. 
Non- parallelism  is  accommodated  within  the  scales  involved. 

The  nonlinear  interaction  is  governed  by  a  partial-differential  system  for  the  vortex 
flow  coupled  with  an  ordinary-differential  equation  for  the  TS-waves’  pressure.  The  solu¬ 
tion  of  this  coupled  system  depends  crucially  upon  so-called  interaction  coefficients  which 
are  functions  of  the  Mach  number;  additionally,  the  TS  waves  are  significantly  affected  by 
the  inclusion  of  compressibility.  It  is  found  that  the  interaction  coefficients,  for  subsonic 
flow,  do  not  differ  significantly  in  nature  from  the  incompressible  ones,  but  as  the  flow 
becomes  supersonic  the  restriction  (for  high  Reynolds  numbers)  that  the  TS  waves  must 
be  directed  outside  the  local  Mach-wave  cone  (Ryzhov,  1984;  Zhuk  &  Ryzhov,  1981)  ex¬ 
cludes  a  particular  flow  solution  which  is  only  possible  for  less  oblique  modes.  The  flow 
properties  point  to  the  second  stages  of  interaction  associated  with  higher  amplitudes. 

It  is  found  that  the  present  formulation  breaks  down  as  the  Mach  number  becomes 
large:  for  then,  even  when  the  presence  of  shock/boundary  layer  interaction  is  neglected, 
the  viscous  sublayers  coalesce  to  form  a  single  boundary-layer.  The  structure  applicable 
in  this  hypersonic  limit  has  been  considered  by  Smith  (1989)  and  Blackaby  (1991). 

The  theoretical  idea  is  basically  that,  if  two  low-amplitude  TS  waves  are  present, 
proportional  to  E\,2  —  exp[i(oX ±  |Z  —  fit)]  say;  then  nonlinear  inertial  effects  produce  the 
combination  E\  =  exp[i/?],=  E^  say,  at  second  order,  among  other  contributions,  i.e. 
a  standing-wave  or  longitudinal-vortex  flow  is  induced.  Here  a,  ft  and  12  are  real-valued 
scaled  wavenumbers  and  frequency,  whilst  X ,  Z  and  t  are  scaled  length-  and  time-scales 
(see  later).  Equally,  the  combination  of  the  vortex  and  one  TS  wave  provokes  the  other 
TS  wave. 

Since  the  Reynolds  number  is  assumed  to  be  large,  the  TS  waves  are  supported  by  the 
triple-deck  structure  (Smith,  1979,89);  however  an  extra  sub-boundary  layer  and  a  further 
streamwise  length-scale  are  necessary  to  capture  their  interaction  with  the  longitudinal 
vortices.  The  present  vortex/ wave  interaction  mechanism  is  very  similar  to  that  of  Hall 
&:  Smith  (1989);  the  difference  is  caused  by  an  error  in  the  latter,  uncovered  by  Smith 
oc  Blennerhassett  (1992).  The  amended  interaction  still  has  the  induced  vortices  lying 
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at  the  top  of  the  lower  deck  but  now  the  forcing  from  the  TS  waves  is  solely  from  an 
inner  boundary  condition.  The  wall-shear  of  the  induced  vortices  modifies  the  wall-shear 
of  the  basic  flow  at  the  same  order  as  the  latter’s  leading-order  nonparallel  correction. 
These  corrections  to  the  wall  shear  force  secondary  TS  waves  in  the  lower-deck,  whilst  the 
amplitude  of  the  primary  TS  waves  here  is  governed  by  an  amplitude  equation  involving 
these  corrections  to  the  wall  shear.  The  behaviour  of  the  primary  TS  quantities  at  the 
top  of  the  lower  deck  then  leads  to  longitudinal-vortex  activity  being  forced  there.  Thus 
the  system  is  truly  interactive:  the  longitudinal  vortices  are  driven  by  the  TS  waves,  the 
amplitude  of  which  is  determined  by  an  amplitude  equation  involving  a  vortex-term. 

We  consider  the  interaction  for  the  case  of  compressible  laminar  flow  over  a  semi¬ 
infinite  plate.  In  the  next  section  the  underlying  boundary-layer  flow  is  outlined  and 
the  triple-deck  structure,  for  such  compressible  flows,  is  reviewed.  In  §3  the  interaction 
equations  are  derived,  and  a  few  special  limiting  cases  of  these  equations  are  considered. 
In  §4  numerical  results  are  reported  and  in  the  last  section  some  conclusions  are  drawn. 

2.  FORMULATION 

We  consider  the  boundary  layer  due  to  high-speed  uniform  flow  of  a  compressible 
over  a  flat  plate.  Suppose  that  L  is  the  distance  from  the  leading  edge,  and  u^,  a^,  p ^ 
and  p,^,  are  the  velocity,  speed  of  sound,  density  and  shear  viscosity  of  the  free  stream 
flow,  then  we  assume  that  the  Reynolds  number,  Re  =  p^u^L/ p^,  is  large.  This  is  not 
unreasonable  as  one  is  already  assuming  the  presence  of  a  boundary  layer.  The  second 
important  parameter  is  the  Mach  number,  Moo  =  which  we  take  to  be  0(1)  for 

the  time  being. 

A  nondimensionalisation  based  on  coordinates  Lx  (where  x  is  in  the  direction  of  flow 
and  y  is  normal  to  plate),  velocities  u^u,  time  Lt/u^,  pressure  p^u^p,  density  p^p, 
temperature  T^T  and  shear  viscosity  p^p  is  adopted,  where  the  subscript  oo  denotes  the 
value  of  the  quantity  in  the  free-stream.  We  assume  that  the  fluid’s  viscosity  and  tem¬ 
perature  are  related  by  Sutherland’s  formula.  Full  details  of  the  Navier-Stokes  equations 
equations  for  compressible  flow;  the  resulting  boundary-layer  equations  and  associated 
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similarity  solutions  can  be  found  in  several  books  and  articles  (eg.  Stewartson,  1964). 
Note  that  only  the  ‘wall-values’  of  the  steady  boundary-layer  flow  solution  occur  within 
scales  considered  in  this  paper;  however,  these  quantities  depend  on  the  choice  of  viscosity- 
temperature  relation  as  well  as  other  factors  such  as  whether  the  plate  is  cooled  and  the 
nature  of  the  external  pressure  gradient. 

2.1  The  3-D  compressible  triple-deck  equations. 

The  underlying  structure,  of  the  vortex-wave  interaction  to  be  considered  later,  is  that 
of  the  three-dimensional,  compressible  Tollmien-Schlichting  (TS)  waves  at  large  values 
of  the  Reynolds  number,  namely  the  three-dimensional  ‘compressible  triple-deck’.  This 
structure  has  been  studied  by,  in  particular,  Zhuk  &  Ryzhov  (1981),  Ryzhov  (1984)  and 
Smith  (1989);  the  two-dimensional,  compressible  triple-deck  theory  was  first  considered 
by  Stewartson  &:  Williams  (1969).  Recently,  Cowley  Sz  Hall  (1990)  and  Duck  &  Hall  (1990) 
have  shown  that  the  theory  can  be  adapted  to  include  the  effects  of  a  shock  for  flow  over 
a  wedge,  and  cylindrical  geometry,  respectively;  whilst  Seddougui,  Bowles  &:  Smith  (1991) 
have  considered  the  effects  of  wall-cooling.  For  definiteness,  we  assume  that  the  flow  is 
supersonic  >1)  during  the  formulation  of  the  interaction  equations;  the  subsonic 

and  other  cases  follow  very  similarly  and  in  §4  results  are  also  presented  for  these  cases. 

In  the  following  scalings,  the  Reynolds  number  is  assumed  to  be  large  whilst  the  other 
factors  are  taken  be  be  0(1).  The  latter  are  introduced  to  normalise  the  resulting  governing 
equations  as  far  as  possible;  however,  the  Mach  number  still  remains  in  the  upper-deck’s 
pressure-disturbance  equation  and  hence  it  appears  in  the  TS-eigenrelation. 

The  streamwise,  spanwise  length-scales  and  the  time-scale,  for  >  1,  are 

3  1  3 

( X  -  Xo  ,Z  -  z0)  =  Re~*K]{X  ,Z),  t  =  Re~l  Xw1,  Hw1  Tw  (M^  -  1  )~*t, 

A'j  =  A^/i^rJ(M^-l)-f;  (2.1a  -  c) 

here  (.ro,  -o )  corresponds  to  the  location  of  the  initial  disturbance  of  the  laminar  base-flow. 
In  the  viscous  sublayer,  or  lower  deck, 

y  =  Re ~  I  A  ttt  ■ -  ?  ft  i  rj  ( Ml  -  1 )  ■ -  5  r, 


4 


(u,v,w)  =  Re-l\ipiT$(ML-l)~i(U,  Rt-?\lnlTjV,  W), 


and 


p  -  Poo  =  Re  *  X^Pw  Tw  J(M^o  -  1)  *P. 

The  resulting  lower  deck  equations  are 

ux  vY  +  wz  =  o, 

+  UUX  +  VUY  +  WUZ  =  -Px  +  Uyy  i 
Wi  +  l/Wx  +  VWy  +  WWz  =  -Pz  +  Wyy , 
Py  =0, 

to  be  solved  subject  to 


(2.2a  -  e) 


(2.3a  -  d) 


U  =  V  =  W  =  0,  on  y  =  0;  t/  -►  (F  +  A(X,  Z,  t)),  as  F  -  oo.  (2.3 e  - /) 


The  main  deck  has 


_  i  I  I 

y  =  Re 


(2.4) 


and  merely  transmits  the  small  displacement  effect,  A ,  across  the  boundary  layer  as  well 
as  smoothing  out  an  induced  spanwise  velocity.  The  displacement,  A,  is  related  to  the 
pressure,  P,  via  a  pressure-displacement  law  stemming  from  matching  solutions  across  the 
three  decks  (see  Smith,  1989). 

The  upper  deck  occurs  where 


y  =  Re~t Ki(Ml  -  l)~h\ 

note  that  the  Mach  number  can  be  scaled  out  of  all  but  the  upper  deck  equations. 


(2.5) 


2.2  Linear  Tollmien-Schlichting  modes. 

The  eigenrelation,  for  linear  supersonic  TS-modes,  is  easily  derived  using  the  triple¬ 
deck  scales  and  equations  discussed  in  the  previous  subsection.  It  can  be  written  in  the 
form 

mV  +  =  MA)«.){4(Afr_1)  -  V.  (2.6a) 
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where  Ai  signifies  t lie  Airy  function, 


k  —  I  A i(q)dq  and  £o  =  (2.6b, c) 

J(o 

Here  o,  fi/2  and  Q  are  the  scaled  wavenumbers  and  frequency  of  the  mode.  The  vortex- 
wave  interaction  to  he  described  concerns  only  neutral  modes,  so  o,  /j  and  Q.  are  all  real. 

In  Figure  1  we  present  the  ‘neutral’  solutions  of  the  eigenrelation  (2.6a),  and  its 
subsonic  counterpart,  for  a  few  (illustrative)  choices  of  the  Mach  number  M^.  Here  (and 
hereinafter)  the  ‘wave  [obliqueness]  angle’  is  defined  by 

0  <  9  =  tan-,(/?/2«)  <  90°.  (2.7a) 


We  see  for  subsonic  values  of  the  Mach  number  (M0 0  <  1)  that  neutral  modes  are  possible 
for  all  wave-angles.  However,  for  increasing  supersonic  Mach  number  values  (M0 o  >  1) 
the  solution  properties  start  to  differ  noticeably,  with  only  an  ever  decreasing  range  of  very 
oblique  TS-wave  propagation  angles,  9 ,  being  possible.  Thus  the  restriction  (2.7a),  which 
can  be  re-written 


6  >  tan-1  (Ml,-  1)2  , 


(2.76) 


is  clearly  evident  in  this  figure.  We  shall  see  later,  once  the  interaction  has  been  formulated 
and  numerical  values  have  been  calculated  for  the  important  interaction  coefficients,  that 
this  restriction  proves  to  be  a  more  significant  ‘compressibility-effect’  on  the  interaction 
than  the  ‘direct’  effect  due  to  the  Mach  number  appearing  in  the  interaction  coefficients. 

Smith  (1989)  gives  a  comprehensive  account  of  the  consequences  of  the  eigenrelation 
(2.6a)  on  the  stability  of  the  flow  to  linear  TS-modes  (note  the  factor  of  2  difference  between 
the  definition  of  ft  here  and  that  used  by  Smith,  1989  and  Blackaby,  1991)  Our  concern 
in  this  paper  is  with  a  vortex-wave  interaction  based  on  these  length-  and  time-scales.  In 
the  next  subsection  we  deduce  the  size  of  additional  x-  and  y-  scales  necessary  to  capture 
this  interaction. 
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2.3  The  interaction  scales. 

In  deriving  the  interaction  scales,  the  same  argument  as  Hall  &  Smith  (1989)  is  fol¬ 
lowed  but  based  on  the  compressible  triple-deck  scales  quoted  in  §2.1.  However,  we  take 
the  coupled  lower-  and  upper-deck  equations  as  our  starting  point,  rather  than  returning 
to  the  compressible  Navier  -Stokes  equations;  this  approach  appears  simpler. 

We  know  that  TS-waves  are  governed  by  the  triple-deck  structure,  and  in  particular 
by  the  unsteady  interactive  boundary-layer  equations  holding  in  the  lower-deck  coupled 
with  the  upper-deck  equations  via  a  pressure-displacement  law.  If  the  three-dimensional 
(3D)  TS-wave  amplitudes  are  comparatively  small,  say  of  order  h  1  relative  to  fully 
nonlinear  sizes,  then  nonlinear  inertial  effects  force  a  vortex  motion  at  relative  order  h2\ 
the  TS-modes  are  taken  to  be  proportional  to 

E\  —  exp[2(«A  4-  ^Z  —  Clt)},  Ei  =  exp[i(oA  —  —Z  —  fit)]  (2.8a,  b) 

and  we  see  that  combinations  yield,  in  particular,  induced  longitudinal-vortex  terms  pro¬ 
portional  to 

£3  =  exp  [i{f)Z)\,  (2.8  c) 

having  only  spanwise  dependence. 

It  can  be  easily  shown  that  spanwise  inertial  effects  (such  as  the  lUW; ,y 1  term  of  the  c  — 
momentum  equation)  decay  slowly  like  \/Y2  resulting  in  the  spanwise  velocity  component 
of  the  induced  vortex  to  grow  logarithmically  like  lnF  (Hall  &:  Smith,  1984,89).  Hall 
&:  Smith  (1989)  introduced  the  concept  of  a  new  sub-layer  (‘the  buffer  layer’)  situated 
within,  and  at  the  top  of,  the  lower-deck,  along  with  a  longer  length-scale  (for  amplitude 
modulation)  to  dampen  down  this  logarithmic  growth.  They  showed  that  the  main  vortex 
activity  was  confined  to  this  region. 

Before  deriving  sizes  for  the  modulation  length-scale  and  the  thickness  of  the  buffer- 

layer,  we  briefly  mention  the  link  between  the  x— scales  present  and  nonparallel  effects.  The 

triple-deck  is  a  local  structure  located  at  nondimensionalised  distance  x  =  xq  from  the 

_  3 

leading  edge.  It  is  short,  its  length  being  0(He  $  K\  )  compared  to  the  0(1)  development 
of  the  underlying  boundary  layer,  and  all  the  A— dependence  of  the  TS  modes  is  taken  to 
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be  in  tlic  E\  and  E>  factors.  The  modulation  of  the  modes  is  assumed  to  he  on  a  longer 
x  -scale,  leaving  the  eigenrelation  (2.Ga)  unaffected.  We  define  this  modulation  x  —  scale, 
.V  say.  by 

x  —  ./■(,  =  +  Rt~  8  A',  A',  Rt~%  Ah  <C  $2  •iC  1  .  (2.9a) 


where  ^2  i 


to  be  determined.  Thus  we  have  multiple  scales  in  x :  formally  we  should  make 


the  replacement 


0  0  R(  -%  I\-  0 

OX  OX  +  i2  &X 


(2.9  b) 


in  the  triple  deck  equations  (2.3). 

At  leading  order  the  wall  shear  \w  =  Au,(x)  is  constant  (with  respect  to  the  A  —  and 
X  —  scales)  but  here  we  wish  to  balance  the  next  order  term  into  the  interaction  equations; 
in  fact,  for  later  convenience  we  have  scaled  the  leading  order  value,  A^xo),  out  of  the 
triple  deck  equations.  As  a  Taylor  expansion  about  the  local  station  x  =  xq. 


Aiu(x)  —  Au^xqXI  +  82X  Xb(xo)  +  •••)> 


(2.19) 


here  At  =  \~l  d\u,/dx  is  0(1)  and  represents  the  first  influence  of  nonparallelism  (stream- 
wise  boundary-layer  growth). 

In  the  buffer  layer,  where  Y  =  iy  say  (6  1),  the  size  of  the  spanwise  velocity 

of  the  induced  vortex  in  the  buffer  layer  is  0(/i2  In  K).  ~  h2  Ini,  leading  to  an  induced 

h2  In 6 ,  by  continuity  (and  noting  that  the  modulation 


1  ■  r  1  ^2  Ae  ^  >  2 

streamwise  velocity  of  order 


—  .  i2JRc*  h2 

is  on  X),  which  alters  the  basic,  shear  Xw  by  a  relative  amount  of  order  — — - —  In  8  and 


this  is  the  same  order  as  the  ‘non-parallel’  Aj-term  if 


Kx  6 


f>2 


62  Re  S’  h2 

A,  T 


Ini  . 


(2.11a) 


The  A— modulation  has  been  introduced  to  damp  the  induced- vortex  velocity  components 
in  the  buffer  layer,  and  this  requires  the  inertial  operator,  Yd/dX,  to  balance  the  viscous 
one,  d2 /dY2,  i.e. 


.3  62  Re  & 


(2.116) 
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One  further  relation  (between  the  unknowns  h,  d  and  6>)  is  i\  1uired  and  results 
from  balancing  the  slower  X  modulation  with  the  nonparallel  effects  too  (i.e.  in  the  x  — 
momentum  equation,  balancing  the  A/,  term  with  Py ),  yielding 

=  Re'&Ki  *  .  (2.11c) 

The  other  two  sizes  now  follow  immediately  from  (2.11a-c): 

i 

6  ~  Rt  IS  A' i  ~  5  and  //  ~  ^i?e~  Ah  ?>  /  In  ^  .  (2.11W.  e) 

The  logarithmic  factor  in  (2.11a,e)  is  important  (Smith  A'  Blennerhass  'tt,  1992):  it  is 
wrong  to  dismiss  it  like  Hall  A  Smith  (19S9).  We  note  that  6  is  large,  whilst  h  and  b2  are 
small,  as  required. 

3.  THE  INTERACTION  EQUATIONS 

The  method  of  deriving  the  interaction  equations  is  identical  to  that  used  by  Hall  A 
Smith  (1989);  essentially,  it  involves  a  standard  weakly-nonlinear  triple-deck  analysis  but 
slightly  complicated  due  to  (i)  the  wall  shear  being  weakly  z- -dependent ,  and  (ii).  the  extra 
(buffer)  layer.  Thus  here  we  present  only  the  briefest  outline  of  the  derivation;  our  main 
concern  in  this  paper  being  the  effect  of  compressibility  on  their  solution  properties.  Fuller 
details  can  be  found  m  Blackaby  (1991)  (see  also  Hall  A  Smith,  1984,89;  and  Smith,  1989). 

3.1  The  equations  and  their  derivation. 

The  interaction  equations  comprise  of  a  set  of  equations  for  the  vortex-terms  in  the 
buffer-layer  coupled  with  an  equation  for  the  TS-wave  pressures  resulting  from  matching 
the  solutions  for  the  waves  between  the  lower-deck  and  the  upper-deck.  The  vortex 
equations  are  forced  by  the  wave  via  a  boundary  condition  whereas  a  vortex  term  explicitly 
enters  the  wave-equation. 

In  the  buffer  layer,  the  vortex  terms,  u 3 3 ( X ,  F ) ,  e3 3 ( A' ,  Y )  and  uyd  A,  F ),  have  the 
following  sizes  relative  to  the  lov.- ei -deck  scales  (2.2) 

{U,  V,  W')vortex  =  h2lnb  {Reib2U^/K\.  <Vf’33,  u>33)E3  +  c.c.;  (3.1) 


9 


here  tin*  notation  of  Hall  A  Smith  (19S9)  has  boon  adhered  to  and  c.c.  denotes  complex 
conjugate.  These  scales,  together  with  the  wall  shear  term  bY  in  tin*  V  expansion,  lead 
to  the  vortex  equations 


m  >  v  ^  "rt.Y  v  —  tdie.-n, 
"\»yT  -  =  °* 

which  must  be  solved  subject  to  the  boundary  conditions: 


(3.2a, b) 


Wy.i(X .  oc )  =  iiriy(X,oo)  =  fi.i:iyy(X  J))  =  0  and  d\{3(A\0)  =  -illK\pu \* ; 

(3.2  c-f) 

here, 

A'  =  1  -  {fi2 /4a2)  (3.3) 

and  pi  i  is  the  amplitude  of  the  TS-waves  which  we  choose  to  be  of  equal  amplitude.  Thus, 
wo  see  that  the  vortex  equations  are  only  forced  by  the  TS-waves  via  a  boundary-condition 
which  matches  the  solution  in  the  buffer  layer  with  that  found  in  the  lower-deck. 

Tlve  desired  equation  for  the  pressure  amplitude  |pn  |  of  the  primary  TS-waves  can  he 
derived  by  solving  the  triple-deck  equations  for  the  primary  and  some  forced  TS-waves; 
after  some  manipulation  we  find  that  it  takes  the  form 


«  \pu(X)\Y  +  b  AfcA|p„(X)j  +  c  fi3 3F(X,0)|p„(X)|  =  0.  (3.4) 


This  equation  was  first  derived  by  Hall  &;  Smith  (1989);  however,  here  the  so-called  com¬ 
patibility  coefficients  a,  b  and  c  are  functions  of  the  Mach  number.  The  presence  of  M 
in  these  coefficients  is  one  of  the  reasons  for  the  solution  properties  for  compressible  flow 
differing  from  those  for  the  incompressible  case.  In  fact,  the  for  supersonic  case 


a  = 


2r,7D£0A 


-2 


3cv 

and  c  =  iD£  0A 


-  iB~? 

-2(2r  17  ,  /?2?,2 


2  B  B 

T+3^  +  1 


6=  - 


2ri7£>£ottA  5 


where 


r  i 


^»(fr>) 

Ai'teoY 


3  +  4 

*(£o) 


„-.B- 

3  47 


5 

"3^ 

(3.5a  —  c) 


r  2  = 


MZoV 


D  =  1  + 


7  =  a2  ,  (P_ 

Ai'(fr)’  7  4’ 
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(3.6a-/) 


g  =  4(A/i  -  1)  ~°2  a',d  A  =  ('«)’■ 

The  coefficients  for  other  flows  can  be  derived  very  similarly  and,  in  fact,  numerical 
results  for  the  subsonic  and  incompressible  cases  are  presented  in  the  next  section.  An 
alternative,  less  physically  motivated,  derivation  of  the  interaction  coefficients  is  outlined 
by  Blackaby  (1991)  who  considers  a  generalised  TS  eigenrelation.  The  quantitative  values 
of  these  coefficients,  and  their  resulting  effect  on  the  interaction  properties  are  considered 
later.  The  (normalised)  interaction  coefficients 


■  c\r  —  Real(b/a)  and  c2r  —  Real(c/a), 


(3.7a,  b) 


are  crucial  to  the  solution  properties;  especially  the  large- X  behaviour. 
It  is  possible  to  derive  a  nonlinear  ‘integro-differential’  equation, 


til  ,  (  \  ~\r  ,  c2r(1'2KAi'{0)  fx  2 

—  +  cirAfcA  +  /  p„  (i/>)(JY 


r  (I)  -47(0) 


J  |pn|2(i/’)('Y  -  0)  Wj|pn|  =  0,  (3.8) 


for  the  pressure  amplitude  |pj  j  |  from  the  previous  equations.  A  similar  equation  has  been 
found  by  Smith  &  Walton  (1989),  in  their  study  of  vortex-wave  interactions. 


3.2  Possible  limiting  forms  for  large-A. 

Let  us  consider  analytically  the  possible  flow  solutions  for  large-.Y.  Hall  &:  Smith 
(1989)  found  four  such  options  for  their  system  of  equations;  however,  the  necessary 
amendments  to  their  work  render  one  of  these  options  is  no  longer  feasible,  namely  that 
of  exponential  growth.  Moreover,  there  is  a  swap  in  the  signature  required  for  the  cru¬ 
cial  quantity  A'c 2r  for  the  finite-distance-blow-up  and  the  algebraic-growth- to-infinity 
eventualities  to  be  possible.  Thus,  the  conclusions,  drawn  later,  for  the  case  of  zero  Mach 
number  are  quite  different  from  those  found  in  Hall  &  Smith  (1989).  In  §4,  numerical 
solutions  of  the  interaction  equations  will  be  presented  and  compared  with  the  large- „Y 
asymptotic  predictions  that  follow. 
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(i)  Option  I:  Finite  distance,  break  up. 


Hall  X  Smith  (19S9)  showed  that  a  possible,  ultimate  behaviour  of  the  nonlinear 
interactive  How,  as  X  increases,  was  that  of  an  algebraic  singularity  arising  at  a  finite 
position,  say  as  X  — *  A’0  ;  this  option  is,  in  fact,  still  possible  for  the  corrected  system 
of  interaction  equations  but  with  some  changes  in  the  details.  The  similarity  forms  they 
proposed  are  appropriate,  apart  from  that  for  the  pressure.  As  A  — >  A"0  ,  the  behaviours 
for  the  interaction  quantities  must,  have  the  forms: 

I/Ml |  ~  (Xo-X )~zp(v),  u-33  ~  (Xo-X r*W(f}), 

1/^,7  ~  (A’u  -  A)_1f(i/),  where  fj  =  Y( A0  -  X )~ * .  (3.9 a  -  d) 

When  these  forms  are  substituted  into  the  interaction  equations,  the  resulting  similarity 
equations  can  be  solved  and  we  deduce  that  we  require  the  quantity 

Kc2r  <  0,  (3.10) 

for  this  option  of  finite-distance  break-up  to  be  a  possible  large-A  state  of  the  vortex- 
wave  interaction.  Note  the  change  of  sign  necessary  for  this  option  to  be  possible;  this 
change  is  due  to  the  modifications  found  necessary  by  Smith  &:  Blennerhassett  (1992). 

The  next  option  that  we  consider  is  less  ‘catastrophic7,  as  far  as  the  laminar  flow  is 
concerned,  with  the  solution  continuing  to  downstream  infinity. 

(ii)  Option  II:  Algebraic  response  at  infinity. 

This  option  is  still  also  possible  with  the  corrected  equations.  The  flow  quantities 
must  have  the  following  forms 

Ipi,i~ u’33~A*4m  u33F ~ xi(ij),  ; f,  =  vx~K  (3.n«-d) 

as  A  — ►  oo . 

It  is  easy  to  show,  from  the  integro-differential  equation  (3.8),  that  we  require 

Kc-ir  >  0  (3.12) 


12 


for  this  option  to  be  possible;  note  that  this  is  also  a  different  result  than  Hall  Smith 
(1989)  found. 

The  third  large- JA  option  proposed  by  Hall  &  Smith  (see  also  Smith  Walton,  1989) 
is  that  of  an  exponential  growth  as  X  — *  oc.  This  option  is  no  longer  possible  as  it  relies 
on  a  forcing  term  in  the  1033-equation  that  is  not  present  in  the  corrected  equations.  A 
further  option,  mentioned  by  Hall  Smith,  that  is  still  feasible;  is  that  of  decoupling  due 
to  linearisation.  Here  the  TS  pressure  disturbance  jpn  |  becomes  very  small,  and  the  vortex 
flow  then  grows  slowly  on  its  own  with  downstream  variable  X  from  its  initial  upstream 
state.  However,  this  option  is  ultimately  unstable  to  the  TS-waves  since  the  nonparallel- 
growth  term,  proportional  to  A&,  will  dominate  the  vortex  skin  friction  ^33  y  (•‘^  ’  0)- 

3.3  The  transonic  and  hypersonic  limits. 


There  are  three  obvious  limiting  cases  to  consider  for  the  value  of  the  Mach  num¬ 
ber;  below  we  consider  the  transonic  and  hypersonic  limits,  —>  1  and  M — < ►  oc 

respectively,  whilst  the  incompressible  case  M 00  =  0  is  considered  in  the  next  section. 

(i)  The  Mach  number  tending  to  unity. 

In  his  study  of  the  eigenrelation  (2.6a),  Smith  (1989)  investigated  various  limiting 
cases,  including  those  of  Moo  — >  1  and  — >  00;  Here  the  ‘transonic  limit’  will  be 

considered;  without  loss  of  generality,  we  suppose  that  the  flow  is  (just)  supersonic  and 
define 

m  =  (M200  -  l)i,  (3.13) 


so  that  m  is  small.  Smith  (1989)  showed  that,  in  this  case,  the  TS-wave  quantities  behave 
like 


(a,/3,f2)~(m  Zot*,rh  ?Q*)  +  -**. 


(3.14) 


Substituting  these  into  the  formula  for  the  interaction  coefficients,  we  find  that 


(cir,c2r)~m  ^(cJr,C2r)  -1 - ,  c*r,c*2r  ~  0(1), 


(3.15) 
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no  that  the  interaction  scales  AM,  the  vortex  disturbances  «33y,  ic33  and  the  TS-  pressure 
amplitude  |/'>n  |  also  need  to  oe  scaled  with  ///: 


A' 


v.uXiV. 


33 


\l>\  1 1)  =  ( >h  «  X  .  in  8  1 


T*  . 


3 

til  8 


33  V 


7 

7//  8  (C 


.13,. 

16  l/'ll 


‘)  + 


(3. 1C) 


These  scales  and  the  resulting  set  of  equations  can  be  used  to  check  numerical  results,  for 
the  general  supersonic  case,  by  providing  a  ‘transonic’  asymptote.  We  do  not  consider  the 
transonic  limit  any  further  here;  Bowles  (1000)  has  considered  transonic  boundary  layer 
transition  but  the  author  is  unaware  of  any  vortex-wave  formulations  for  the  transonic 
regime. 

(ii)  The  large  Mach  number  limit. 

Another  limiting  case  that  Smith  (1080)  went  on  to  investigate  was  the  so-called 
hypersonic  limit  when  M0 0  1;  this  limiting  case  leads  to  some  interesting  consequences 
for  the  whole  triple-deck  structure.  Thus  it  would  be  most  instructive  to  consider  the 
same  limit  here,  as  our  interaction  structure  is,  of  course,  dependent  to  a  very  great  extent 
on  the  underlying  triple-deck  scales. 

First,  let  us  recap  those  results  of  Smith  (1080)  which  are  relevant  here,  before  going  on 
to  investigate  the  result  of  increasingly  large  Mach  number  on  the  interaction-coefficients, 
equations  and  length-scales.  For  M0 0  1  the  main  features  revolve  around  the  small 

regime 

(a,  /?,  fl)  ~  (MZ? MZ?  V,  M~l  ft)  +  •  ■  • ,  (3.17) 

where  d,  0  and  ft  are  0(1).  Since  a,0  and  ft  appear  in  the  interaction  coefficients  (3.7), 
it  is  necessary  to  rescale  the  coefficients  as  follows: 

_  3 

(cir,c2r)  =  Moo2(cir,c2r)  +  •••,  where  C],c2~0(l).  (3.18a, b) 


It  also  is  necessary  to  rescale  the  quantities  appearing  in  the  interaction  equations, 

(X,F.u„7,«.„,|p„|)  =  (Ml/'X.M'J'Y.MU'inr.MU'An,  Af-V8 1£„|). 

(3.18c- <7 ) 
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A  set  of  interaction  equations  for  this  hypersonic-limiting  case  can  be  derived  (Black- 
aby,  1991).  These  appear  exactly  the  same  as  the  general  interaction  equations  which 
obscures  the  fact  that  the  whole  multi-layered  boundary-layer  structure  is  radically  al¬ 
tered  as  the  Mach  number  increases.  It  was  shown  by  Blackaby  (1991),  based  on  the  theory 

i  _  3 

of  Smith  (1989),  that  as  Moo  Z  the  triple-deck  streamwise  length-scale,  Rt  % I\\X , 
rises  to  become  0(1)  in  size;  implying  that  a  normal-mode  decomposition  is  no  longer 
rational  because  the  nonparallelism  of  the  underlying,  growing  boundary-layer  is  now  a 
leading-order  effect.  Further,  in  this  limit  it  was  shown  that  the  lower-deck  thickens  to 
coalesce  with  the  main-deck. 

This  collapse  of  the  underlying  compressible-triple-deck  structure,  as  the  Mach  num¬ 
ber  increases,  will  obviously  occur  for  the  large  Mach  number  behaviour  of  the  vortex-wave 
interaction  being  considered.  However  the  large  Mach  number  destiny  of  buffer-layer  (in 
particular,  its  thickness)  and  the  amplitude-modulation  scale  remain  to  be  established. 
Intuitively,  as  the  buffer-region  is  ‘sandwiched’  between  the  lower-  and  main-decks  which 
merge  into  a  single  viscous  layer  in  this  limit,  we  would  also  expect  the  buffer- region  to 
collapse  into  the  same  viscous  layer.  Similarly,  as  the  modulation-scale  is  ‘sandwiched’ 
between  the  triple-deck’s  streamwise  length-scale  (which  emerges  as  0(1)  in  this  Mach 
number  limit)  and  the  0(l)-length-scale  of  the  underlying  flow,  we  would  expect  that  the 
modulation-scale  also  lengthens  to  that  of  the  underlying  base  flow  (as  Moo  Z  Re^).  We 
now  show  that  these  suspicions  are  correct,  by  formally  considering  the  large  Mach  number 
properties  of  the  scales  involved. 

Recall  that,  in  the  streamwise  direction,  we  have  the  multiple  scales, 

dx  — >  dx  +  "  ReZK^'dx; 

necessary  to  capture  the  vortex-wave  interaction.  The  quantities  A'i  and  6?  are  as  defined 
by  (2.1c)  and  (2.11c),  respectively.  In  the  large  Mach  number  limit,  we  have  seen  that 

_  3  _  3 

~  Moo4  whilst  dx  ~  or  ~  Mx°  , 
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so  that  the  unsealed  length  scales,  Lw  and  Lv  say,  of  the  TS  waves  and  the  modulation 
of  the  induced  vortices,  respectively,  are 

3  3  3 

Lw  ~  Re  8  7\  i  <C  1  and  Lv  ~  b2 M<J,  1. 


The  Sutherland  temperature-viscosity  relation  leads  to. 


Ad 


15 

Ml. 


(3.19a) 


and 


so 


3  27  | 

Lw  ~  Re~  »  ,  /  0{  1),  as  f  Re  ? 


As  far  as  the  amplitude-modulation  scale  is  concerned,  we  find  that 


3  13  2  T 

Lv  ~  S2M^  ~  Re~t$K^M^i  ~  Re~  T6  S  0(1),  as  Moo  Ac 5 


(3.196) 


Thus,  as  predicted  earlier,  this  modulation  scale  does  indeed  rise  to  0(l)-size  in  this  limit 
of  the  Mach  number. 

Now,  let  us  consider  the  buffer-region;  it  lies  at  the  top  of  the  lower-deck,  where 
the  lower-deck  normal- variable  Y  =  6Y  and  6  is  defined  by  (2. lid).  For  large  Mach 
numbers,  we  have  found,  (3.18d),  that  the  buffer-region  is  characterised  by  the  location 

_  l _ 

where  Y  =  MS*  Y  ~  0(1).  Thu  -  the  buffer-region  lies  where 
Y  ~  SmI  ~  ReTZK~K\li  ~  ReTEMn™ 

->  O(ReTS)  ~  0(m1),  as  /  Re? .  (3.20) 

Note  that  for  large  Mach  number,  the  lower-deck  variable,  Y ,  also  scales  on  M^;  in  fact 
1 

Y  ~  M£o  —  hence  from  this  and  (3.20)  we  deduce  that  the  buffer-layer  merges  with  the 
lower-deck,  which  in  turn  coalesces  with  the  main-deck.  Thus  the  three  sub-boundary- 
layers,  present  for  M ^  ~  0(1),  have  all  merged  into  one  single  viscous  layer. 

Summarising,  when  Moo  — >  Re?  the  four-layered,  short-scaled  structure  underlying 
the  vortex-wave  interaction  collapses  into  the  two-tiered,  long  structure  found  by  Smith 
(1989)  and  considered  by  Blackaby  (1991). 
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4  RESULTS  AND  DISCUSSION. 

This  study  was  motivated  by  the  desire  to  find  out  what  changes  to  the  theory,  pre¬ 
dictions  and  conclusions  of  the  original  work  by  Hall  &  Smith,  are  brought  about  by  the 
inclusion  of  compressibility  effects.  However,  ironically,  the  changes  brought  about  by  the 
correction  of  the  former  turn  out  to  be  more  significant.  For  this  reason,  and  for  later 
comparison,  the  new  results  for  incompressible  flow  will  also  be  presented  in  §4.2.  In  the 
following  subsection,  we  show  how  the  interaction  equations  can  be  ‘normalised’  so  that 
their  solution  depends  merely  on  initial  conditions  imposed  and  the  sign  of  A'c2r.  Numer¬ 
ical  solutions  of  the  normalised  interaction  equations  are  presented  for  both  choices  of  the 
sign  of  A'c2r. 

4.1  The  interaction  equations  renormalised. 

In  §3.2  we  considered  possible  limiting-forms,  for  solutions  to  the  interaction  equa¬ 
tions  as  A  — »  oo,  and  found  that  the  sign  of  the  quantity  Kc2r  was  crucial  in  deciding 
whether  particular  limiting  forms  were,  in  fact,  possible.  This  suggests  that  the  interaction 
equations,  (3. 2, 3. 4),  can  be  renormalised.  This  being  desirable,  we  investigated  further  and 
found  this  was,  indeed,  the  case. 

Writing 

X  =  |clrA6|-7X%  F  =  |clrA6|-£r*,  w33  =  -iPK\c3\ ~'W\ 

IPnMcl-ip*  and  “33y  =  -p2K\c\Ab\~^\cz\~'T*,  (4.1a  -  e) 

where 

c3  = -A'c2r/?2/|c,rAjj§,  (4.2) 

leads  to  the  normalised  system 

Wy.y.  -  Y*W*X.  =  0,  Ty.y.  -  Y*  T*X.  =  W* 
and  P*x.  +  [sgn(clrA6)A*  -  sgn(A'c2r)r*(A*,  0)]  P*  =  0,  (4.3a  -  c) 


17 


which  must  be  solved  subject  to  initial  conditions  (at  X*  =  0),  together  with  the  boundary 
conditions 


W*(A*,  oo)  =  r*(A*,  oo)  =  ry-.(A'*,  0)  =  0  and  VT(X\  0)  =  P*2(A* ).  (4.3d  -  g) 


Thus,  the  interaction  equations  (and  hence  their  solutions)  are  dependent  only  on  the 
initial  conditions  imposed;  sgn(cirAj)  and  sgn(AY2r).  In  all  the  numerical  calculations 
carried  out,  it  was  found  that  cir  >  0,  whilst  <  0  for  a  growing  ‘Blasius-similarity- 
variable-type’  boundary  layer  —  appropriate  to  the  present  study,  if  we  assume  that  there 
is  no  significant  wall-cooling  or  pressure-gradient  effects.  We  therefore  set  sgn(cirA6)  =  —  1 
and,  apart  from  consistent  initial  conditions,  the  only  parameter  remaining  is  sgn(A'c2r)- 
Thus,  with  hindsight,  it  is  not  surprising  that  the  (predicted)  solution  properties  for,  large- 
A,  depend  crucially  on  the  value  of  sgn(A'c2r)-  Recall  that  earlier,  in  §3.2,  we  noted  the 
following  predictions: 


sgn(AY2 


-2r)  | 


>  0  :  Algebraic  response,  as  X  — ►  oo 

<  0  :  Finite-distance  break-up,  as  X 


X 


(4.4) 


o  > 


for  the  behaviour  of  the  solutions  to  the  interaction  equations. 

To  check  these  predictions,  the  normalised  system,  (4.3),  was  solved  numerically;  for 
both  possible  values  of  sgn(A'c2r),  and  for  different  (consistent)  initial  conditions.  The 
large- A  (A*  1)  properties  of  the  solutions  were  found  to  depend  solely  on  sgn(A'c2r); 

the  initial  conditions  were  found  to  affect  only  the  initial  development  of  the  imposed 
disturbances.  The  equations  were  solved  by  taking  ‘central  differences’  in  Y*  and  ‘forward 
differences’  in  A*  (following  the  method  of  Hall  &  Smith,  1989);  the  appropriate  numerical 
checks  were  performed. 

In  Figures  2a, b,  we  present  typical  results  for  both  values  of  sgn(Ac2r)-  In  both  of 
these  computations  the  system  was  initialised  at  A*  =  —  1  (upstream  of  the  neutral  TS 
point)  using 


P*  =  p0*,  W *  =  P0*2(  1  +  y*2)exp[-r*2],  T*  =  (l-^-  +  r*2)exp[-r*2],  (4.5a  -  c) 


’■*2  • 


'*2 


Pn*2 


*2' 


'*2l 


with  Pq  =  0.1.  Note  that  this  initial  state,  which  is  consistent  with  the  interaction 
equations  plus  boundary  conditions,  corresponds  to  a  ‘mixed’  wave/vortex  state.  Moreover, 
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we  set*  from  the  ‘forcing’  boundary  condition  {t*g.  (4.3g))  that  admissible  initial  states 
cannot  consist  of  just  TS  waves  alone;  the  longitudinal  vortices  must  initially  be  present. 
It  appears  to  the  author  that  the  initial  states  used  by  Hall  Smith  (1989)  (see  their  section 
5;  particularly  figures  2  5)  are  inconsistent  with  their  system  of  interaction-equations  plus 
boundary  conditions;  they  do  not  appear  to  satisfy  the  boundary  conditions.  In  their 
study  of  vortex/wave  interactions,  Smith  Walton  (1989)  do  not  comment  on  the  initial 
conditions  they  choose. 

Returning  to  Figures  2a, b,  we  see  that  these  numerical  results  are  in  full  agreement 
with  the  theoretical  large- X*  predictions,  (4.4).  Thus,  in  the  following  subsections,  it 
is  sufficient  to  merely  calculate  values  of  sgn(A”c2r)  in  order  to  determine  the  solution 
properties  for  large-.Y*;  these  being  of  principal  interest. 

4.2  The  incompressible  case  (M0 0  =  0). 

In  their  study,  Hall  &  Smith  (1989)  considered  ‘this’  vortex-wave  interaction  for  in¬ 
compressible  boundary-layer  flow.  They  cleverly  deduced  the  scales  and  formulated  the 
interaction;  unfortunately,  they  made  two  unrelated  errors  in  their  analysis,  both  of  which 
have  a  significant  effect  on  the  results  and  conclusions.  The  first  of  these  errors  concerning 
the  missing  logarithmic  term  in  the  interaction  scales  (most  kindly  pointed  out  to  the 
current  author  by  Dr.  P.  Blennerhassett  and  Prof.  F.T.  Smith),  leads  to  a  simpler  system 
of  interaction-equations,  as  well  as  leading  to  changes  in  the  possible  large-X  states  and 
the  necessary  parameter  values  for  them  to  be  possible.  The  second,  the  term  ue^A $zY 
missing  from  the  left-hand  side  of  Hall  &  Smith’s  (1989)  equation  (3.9b),  was  spotted  by 
the  current  author  and  leads  to  a  corrected  form  for  c,  and  hence,  a  corrected  value  for 
the  crucial  quantity  C2r- 

The  interaction  coefficients  «,  6,  are  as  given  by  Hall  &  Smith  (1989);  whilst  the 
corrected  form  for  c  was  given  by  Blackaby  (1991)  (see  also  the  Appendix  A  of  Smith  & 
Blennerhassett,  1992).  In  Figure  3  the  new  numerical  values  for  the  important  interaction 
quantities,  nr  and  c2r  are  plotted,  versus  TS-wave  obliqueness  angle  0\  recall  that,  for 
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incompressible  flow,  all  such  wave  angles  are  possible.  Note  that  cir  >  0  for  all  0;  whilst 
C2r  has  one  zero,  at  0  ~  32.21°.  Recalling  the  definition  of  A', 


K 


=  1  —  (/32/4a2)  —  l  —  tan2 


>  0 
<  0 


:  if  0  <  45° 

:  otherwise, 


we  see  that,  when  M ^  =  0, 


:  if  32.21°  <  0  <  45° 
:  otherwise. 


Thus,  from  this  last  result  and  the  numerical  calculations  described  in  §4.1,  we  deduce 
the  following:  (i)  if  32.21°  <  0  <  45°  then  the  solution  to  the  interaction  equations  will 
‘blow-up’  in  a  finite-distance;  otherwise  (ii)  the  solutions  will  grow  slowly  (far  slower  than 
the  linear  TS-solutions  if  there  were  no  vortices  present),  with  amplitudes  proportional  to 
algebraic  powers  of  X* ,  as  X*  — >  oc.  Note  that  these  conclusions  are  quite  different  from 
those  of  Hall  &:  Smith  (who  concluded  that  the  ‘finite-distance  break-up’  option  was  most 
likely,  apart  from  the  small  range  45°  <  0  <~  50°  where  an  ‘exponential-growth’  option 
was  favoured).  Thus,  the  theoretically-exciting  ‘finite-distance  break-up’  option  is  now 
the  exception,  rather  than  the  rule. 

4.3  The  subsonic,  supersonic  and  hypersonic  cases. 

For  subsonic  (M<x>  <  1)  and  some  supersonic  (1  <  Moo  <~  1.15)  flows,  the  properties 
of  the  interaction-coefficients  were  remarkably  similar  to  those  found  for  the  incompress¬ 
ible  case  i.e.  graphs  of  Cir,C2r  against  0  appear  very  similar  to  Figure  3.  However,  the 
TS-wave  angle  restriction  (2.7b)  is  found  to  have  a  far  more  significant  effect  for  ‘more’ 
supersonic  flows  —  essentially  it  can  be  regarded  as  preventing  wave-angles  that  would 
allow  sgn(A'c2r)  <  0,  corresponding  to  the  finite-distance  break-up  option.  This  is  illus¬ 
trated  more  clearly  in  Figure  4  where  the  results  are  summarised;  we  see  that  the  0  -•  Moo 
plane  splits  into  four  regions  (labelled  I  -  IV,  as  shown).  Region  IV  corresponds  to  the 
‘barred’  area,  where  no  neutral  TS-modes  are  possible.  We  see  how  the  border  of  this 
region  acts  as  an  ‘abrupt  cut-off’  to  the  larger-Moo  extent  of  Region  II  (  finite-distance 
break-up  option).  This  is  so  much  so  that,  for  Mach  numbers  above  \/2,  the  possibility 
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of  finite  distance  break  up  has  gone.  Thus  summarising,  in  the  subsume  ease  the  results 
are  almost  identical  to  the  incompressible  case:  whereas,  in  general,  the  finite  distance 
break  up  eventuality  is  not  possible  for  supersonic  Hows,  mainly  due  to  the  severe  cut  oH 
restriction  for  large  Reynolds  numbers.  To  illustrate  the  last  point,  in  Figure  o  we  have 
plotted  <•■),■  versus  ,1  and  8  for  M =  3  note  (i)  that  then*  is  no  zero  for  and  (ii). 
the  very  oblique  wave  angles  encountered  (so  that  A  is  always  negative  and.  hence,  /vis, 
is  always  positive). 

The  last  set  of  results  that  we  present  are  for  hypersonic  How  over  a  wedge,  as  consid¬ 
ered  by  Cowley  &:  Hall  (1990),  in  which  a  shock  is  fitted  into  the  upper- -deck  (at  y  =  yH. 
where  y  is  the  normal- variable  of  the  upper  deck),  leading  to  a  modified  form  of  Smith's 
hypersonic  TS-eigenrelation  (the  reader  is  referred  to  the  paper  by  Cowley  <k:  Hall  for  all 
details  of  the  formulation).  In  Figure  G.  we  present  results  for  the  first  (lowest)  neutral 
curve  for  the  case  y,  =  1;  here  Of//,  ficM  ~  0(1)  correspond  to  atjj  in  the  notation  of 
that  paper.  It  is  sufficient  to  note  that,  in  our  notation, 


£ 
2  a 


Moo 


0CH 

C*CH 


>  1, 


and  so  the  waves  they  consider  are,  in  general,  very  oblique.  Of  particular  interest  here  is 
the  (small)  interval  where  C2r  >  0,  so  that  Kc'j.r  <  0,  corresponding  to  the  finite-distance 
break-up  option;  this  is  an  effect  of  the  shock.  No  such  interval  is  found  for  the  ‘higher’ 
neutral  curves;  this  interval  appears  to  be  a  feature  of  the  ‘lowest’  neutral  curves  only  (for 
each  choice  of  ys )  and  corresponds  to  ‘crossing’  the  ‘divide’  c/ch  —  0c  H  ■ 

Finally,  we  report  that  for  the  ‘hypersonic  and  transonic’  limiting  cases  mentioned  in 
§3.3,  the  numerical  results  and  the  predicted  asymptotic  behaviours  (for  the  interaction 
coefficients)  were  in  extremely  good  agreement. 


5  CONCLUSIONS 


Many  of  the  conclusions  of  Hall  &  Smith  (1989)  carry  over  to  the  present  study 
and  so  we  concentrate  on  compressibility-related  aspects  here.  In  this  paper  it  has 
been  demonstrated  that,  within  the  triple-deck  framework  (Re  1),  pairs  of  small- 
amplitude  Tollmien-Schlichting  waves  and  longitudinal  vortices  can  interact,  leading  to 
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mutual  growth.  We  have  seen  that  two  possible  'eventualities’,  for  the  downstream  evolu¬ 
tion  of  the  interaction,  exist;  one  in  which  the  solutions  grow  relatively  slowly  as  A  — +  dc: 
whilst  the  other  terminates  at  a  finite-distance  in  a  ‘break-up'.  Further,  we  have  seen  that 
the  latter  is  no  longer  possible,  in  general,  for  supersonic  flows  (Rt  1). 

The  interaction  has  been  considered  for  all  ranges  of  the  Mach  number;  corrected 
results  for  the  incompressible  case  have  been  presented  and  the  main  effect  of  compress¬ 
ibility  is  through  its  impact  on  the  TS  waves  via  their  governing  mechanism,  the  triple  deck 
structure.  In  the  transonic  and  hypersonic  limiting-cases  the  interaction  modulation  scale 
A'  must  be  rescaled;  in  the  transonic  limit  this  modulation  scale  shortens,  whilst  in  the 
hypersonic  limit  the  opposite  is  true.  The  investigation  of  such  vortex/ wave  interactions  in 
transonic  and  hypersonic  flows  (not  their  ‘limits’)  should  prove  interesting;  note  that  the 
former  flow  has  been  studied  by  Bowles  (1990),  whereas  the  latter  flow  regime  has  been 
considered  by  Blackaby  (1991). 

Other  effects  which  could  be  incorporated  into  the  present  theory  include  pressure- 
gradient  effects;  wall  cooling  effects  (see  Seddougui.  Bowles  k  Smith,  1991);  cylindrical 
geometry  (see  Duck  k  Hall,  1990)  and  spamvise-modulation  (cf  Smith  k  Walton.  1989). 
Finally,  we  note  that  Hall  k  Smith  (1991)  and  Walton  k  Smith  (1992)  consider  the  pro¬ 
portion  of  ‘strongly  nonlinear’  TS-wave/vortex  interactions  corresponding  to  larger  wave 
amplitudes  than  those  considered  in  this  paper. 
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Figure  2b. 


Figure  2b.  Numerical  solution  of  the  interaction  equations  (4.3)  with  A'c2r  >  0. 


Figure  4. 
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